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The investigation of the coupled atmosphere-ocean system is not only sci- 
entifically challenging but also practically important. 

We consider a coupled atmosphere-ocean model, which involves hydrody- 
namics, thermodynamics, and random atmospheric dynamics due to short time 
infiuences at the air-sea interface. We reformulate this model as a random dy- 
namical system. First, we have shown that the asymptotic dynamics of the 
coupled atmosphere-ocean model is described by a random climatic attrac- 
tor. Second, we have estimated the atmospheric temperature evolution under 
oceanic feedback, in terms of the freshwater fiux, heat flux and the external 
fluctuation at the air-sea interface, as well as the earth's longwave radiation 
coefficient and the shortwave solar radiation profile. Third, we have demon- 
strated that this system has finite degree of freedom by presenting a finite 
set of determining functionals in probability. Finally, we have proved that the 
coupled atmosphere-ocean model is ergodic under suitable conditions for phys- 
ical parameters and randomness, and thus for any observable of the coupled 
atmosphere-ocean flows, its time average approximates the statistical ensemble 
average, as long as the time interval is sufficiently long. 



1. Geophysical background 

The coupled atmosphere-ocean system defines the environment we live. Random- 
ness or uncertainty is ubiquitous in this coupled, complex, multiscale system: For 
example, stochastic forcing (wind stress, heat flux and freshwater flux), uncertain 
parameters, random sources or inputs, and random boundary conditions. 
Mathematical models are a key component of our understanding of climate and 
geophysical systems. It is our belief that the fidelity of these models to nature can 
greatly benefit through the inclusion of stochastic effects. Taking stochastic effects 
into account is of central importance for the development of mathematical models 
of many phenomena in geophysical and climate flows. 

We consider a two-dimensional coupled atmosphere-ocean model in the latitude- 
depth plane, with atmospheric dynamics highly simplified, i.e., the atmospheric 
dynamics is described by an energy balance model. The oceanic dynamics is de- 
scribed by the Navier-Stokes equation in vorticity form and the transport equations 
for heat and salinity. The energy balance model is under random impact due to, 
for example, eddy transport fluctuation, stormy bursts of latent heat, and flickering 
cloudiness variables. So this coupled atmosphere-ocean model consists of stochastic 
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and deterministic partial differential equations, together with air-sea flux or Neu- 
mann boundary conditions. We will reformulate this model as a random dynamical 
system. 

The ocean and the atmosphere are constantly interacting through the air-sea ex- 
change process. The ocean moves much slower than the atmosphere does. It is 
generally believed that the ocean plays an important role in the global climate 
dynamics in relatively long time scales, due to ocean's large capacity of holding 
and transporting huge amount of heat or cold around the globe |27| . However, a 
complete quantitative understanding or estimate for ocean's impact on climate is 
lacking. A particular issue is: How does the ocean affect or provide feedback to 
the air temperature, which is the most important climate quantity we are usually 
concerned about? This is a direct impact of the ocean on the climate. It is de- 
sirable to predict or estimate this feedback in the context of our simple coupled 
atmosphere-ocean model. 

The existence and interpretation of climatic attractors have been controversial and 
have caused a lot of debate JHI ■ A low dimensional climatic attractor was regarded 
as an indication that the main feature of long-time climatic evolution may be viewed 
as the manifestation of a deterministic dynamics. We will consider random climate 
attractors, and the long time regimes that such attractors represent still carry the 
stochastic information of the geophysical flow system. We will also investigate 
the finite dimensionality of the asymptotic dynamics by checking the determining 
Junctionals in probability. 

In a special case of physical parameters and random noise, we obtain a random at- 
tractor which is defined by a single random variable. This random variable attracts 
all other motions exponentially fast. This random variable corresponds to a unique 
invariant measure, which is the expectation of the Dirac measure with the random 
variable as the random mass point; see [21 ■ In this case, the coupled atmosphere- 
ocean model is ergodic, and thus for any observable of the coupled atmosphere-ocean 
flows, its time average approximates the statistical ensemble average, as long as the 
time interval is sufficiently long. 

In the next section, we present the coupled atmosphere-ocean model, and discuss 
the well-posedness of this coupled model in §3. Then we investigate the dissipa- 
tivity property in §4. This property is the basis of the asymptotic behavior of the 
coupled system to be considered in §5: atmospheric temperature evolution (with 
oceanic feedback), random attractors, finite dimensionality and ergodicity. Finally, 
we summarize these results in §6. 



2. A COUPLED ATMOSPHERE-OCEAN MODEL 

We consider a zonally averaged, coupled atmosphere-ocean model on the merid- 
ional, latitude-depth {y, 2;)-plane as used by various authors [23 EH El El El • It is 
composed of a one-dimensional stochastic energy balance model proposed by North 
and Cahalan |2()| . for the latitudinal atmosphere surface temperature Q{y,t) on 
domain < y < 1, together with the Boussinesq equations for ocean dynamics in 
terms of vorticity q{y, z, i), and transport equations for the oceanic salinity S{y, z, t) 
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and the oceanic temperature T{y, z, t) on the domain D = {(y, z) : < y, z < 1}: 
e* ^Oyy - (a + 9) + Saiy) - b{y){Soiy) + © " T{y, 1)) + w, 
qt + J{q, ^) =PrAg + Pr ' ^&{dyT - dySy), 

St + J{S, ij) ^AS, 

where 

q{y,z,t) = -All) 

is the vorticity, a is a positive constant parameterizing the effect of the earth's 
longwave radiative cooling, Sa{y) and So{y) are empirical functions representing 
the latitudinal dependence of the shortwave solar radiation, b{y) is the latitudinal 
fraction of the earth covered by the ocean basin, Pr is the Prandtl number and Ra 
is the Rayleigh number. The first equation is the energy balance model proposed 
by North and Cahalan [20]. The fluctuating forcing w{y^t) may arise from, for 
example, eddy transport fluctuation, stormy bursts of latent heat, and flickering 
cloudiness variables. This forcing term is usually of a shorter time scale than the 
response time scale of the large scale oceanic thermohaline circulation. So we neglect 
the autocorrelation time of this fluctuating process as in [201 • We thus assume the 
noise is white in time. The random white-in-time noise w{y^t) is described as 
the generalized time derivative of a Wiener process w{y,t) with mean zero and 
covariance operator Q. Moreover, J{g, h) — gxhy — gyh^ is the Jacobian operator 
and A — dyy -I- dzz is the Laplacian operator. All these equations are in non- 
dimensionalized forms. 

Note that the Laplacian operator A in the temperature and salinity transport 
equations is presumably dyy + ^5"^ dzz with 5 being the aspect ratio, and khi 
the horizontal and vertical diffusivities of heat/salt, respectively. However, our 
energy-type estimates and the results below will not be essentially affected by taking 
a homogenized Laplacian operator A = dyy + dzz- All our results would be true 
for this modified Laplacian. The effect of the rotation is parameterized in the 
magnitude of the viscosity and diffusivity terms as discussed in [31] ■ 
The no-flux boundary condition is taken for the atmosphere temperature Q{y,t) 

ey{o,t) = eyii,t) = o. 

The fluid boundary condition is no normal flow and free-slip on the whole boundary 

^^0, g = 0. 

The flux boundary conditions are assumed for the ocean temperature T and salinity 
S. 

At top z = 1, the fluxes are specified as: 

(2) dzT{y,l) = So{y) + e{y)-T{y,l), dzS{y,l) ^ F{y), 

with F{y) being the given freshwater flux. 
At bottom z — 0: 

dzT^dzS = 0. 
On the lateral boundary y S {0, 1}: 

dyT = dyS = 0. 

The stochastic partial differential equation for air temperature O in Q is only 
defined on the air-sea interface (0 < ?/ < 1) and it may be regarded as a dynamical 
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boundary condition. The boundary condition Q involves a coupling between the 
atmospheric and oceanic temperature at the air-sea interface. 
The deterministic version of this model was studied in |14j . Now we look at the 
well-posedness of this coupled atmosphere-ocean model and then investigate its 
random dynamics. 

3. Well-posedness 

In this section we will show that defines a well posed model. In particular, we 
can show that has a unique solution. Without such a property it would not be 
possible to make predictions from the model numerical simulations or investigate 
the stability behavior. 

Now we are going to re-formulate the model such that appropriate tools of the the- 
ory of random dynamical system can be applied to analyse the coupled atmosphere- 
ocean model under a random wind forcing. For the following we need some tools 
from the theory of partial differential equations. 

Let W2{D) be the Sobolev space of functions on D with first generalized derivative 
in L2{D), the function space of square integrable functions on D with norm and 
inner product 

\W\\l2 — (^J \u{x)\^dD^ , {u,v)l2— J u{x)v{x)dD, u, v e L2{D). 
The space W2{D) is equipped with the norm 

IWWwi = + \\dyU\\L2 + \\dzU\\L2- 

Q 

Motivated by the zero-boundary conditions of q we also introduce the space wl {D) 
which contains roughly speaking functions which are zero on the boundary dD of 
D. This space can be equipped with the norm 

(3) \W\\wl{D) = \\dyU\\L2 + \\dzU\\L2- 

Similarly, we can define function spaces on the interval (0, 1) denoted by £2(0, 1) 
and 1). 

Another Sobolev space is given by Wi{D) which is a subspace of W2{D) consist- 
ing of functions u such that udD = 0. A norm equivalent to the 11^2^ -norm on 
Wi{D) is given by the right hand side of (Q. For functions in L2{D) having this 
property we will write L2{D). 

Since D has a Lipschitz continuous boundary dD there exists a continuous trace 
operator. 

j9D:W^{D)^H^dD). 

Here H^{dD) is a boundary space, see Adams ^ or below. Similarly, we can 
introduce trace operators that map onto a part of the boundary of dD for instance 
for the subset {{y,z) e D\z — 1} oi D. For this mapping we will write 

(4) ^.^i:WI{D)^h"^{0,1). 
The adjoint operator 

is also continuous. Note that ' denotes the dual space for a given Banach space. 
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Our intention is now to formulate the problem |^ with the non-homogeneous 
boundary conditions in a weak form. For convenience, we introduce the vector 
notation for unknown geophysical quantities 

(5) u^ie,q,T,S). 



We now take the linear differential operator from 



Au = 



( -dlyQ + (1 + b{v))Q \ 
-PrAg 
-AT 

V / 



Remind that the function 1 > h{y) > 0. A is defined on functions that are suf- 
ficiently smooth. We also have the following boundary conditions from Section 

m 

dyQiQ) = 9,6(1) - 0, 
^\aD = (^, 9|ai3 = 0, 

d,T{v, 1) = So[v) + e(y) - T{y, 1), d,S{y, 1) - F{y), 
d,T{y,0)^d,S{y,0)^0, 

dyT{Q, z) = dyS{Q, z) = 0, 9yT(l, z) = d^S{l, z) = 0. 
We will assume that So, Sa and F e ^2(0, 1). Note that 

— j Sdydz = j F(y)dy = constant. 

It is reasonable (see ^) to assume that 



(6) 



F{y)dy = 0, 



and thus Sdydz is constant in time and we may assume it is zero: 



Sdydz = 0. 



D 



Thus we have the usual Poincare inequality for S. Unfortunately, this is not the 
situation for T. However, we can derive the following Poincare inequality 



(7) 



|Tf <2||7,=ir||i. + 4||vrf. 



as in Temam [301, Page 52. 

We introduce the phase space for our geophysical quantities H = ^2(0, 1) x L2{D) x 
L2{D) X L2{D) with the usual L2 inner product and V = ^^(0,1) x W\{D) x 
W2{D) X W2{D). For another sufBciently smooth functions v = {Q,q,f,S), we 
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can calculate via integration by parts 
(8) 



{Au, v)h 



dyQdyedy+ / (1 + 6)69^2/ 



y^^y 
(0,1) "'O 

-Pr / Vq-VqdD 

JD 



VT-VfdD- / {So{y) + Q{v)-lz=iT{v,z)){^,^if{y,z))dy 



D 

VS-VSdD- f' F{y){j,=,S{y,z))dy. 

ID JO 

Hence on the space V we can introduce a bilinear form a(-, •) which is continuous, 
symmetric and positive 

d{u, v)^ [ dyQ ■ dyQdy + / (1 + b)QQdy + Pr [ Vq- VqdD 

J(0,l) JO JD 



L 



VT ■ \/TdD + Co / jz^iTj.^iTdy + VS ■ VSdD 

D Jo JD 

for some cq > 0. The other terms from © will be considered separately. This 
bilinear form defines a unique linear continuous operator A : V V such that 
{Au,v) — d{u,v). We can see that the bilinear form d{u,v) is positive, using the 
Poincare inequality (0 . According to , we now introduce the nonlinear operator 

F{u) := Fi{u) + F2[u) where 

( \ 

-J{q,^) 

-J{T,^) 



Fi{u)[y,z\ = 



and 



F2{u)[y,z] 



J 



( ^a + Sa{y)^b{y){So{y)-i.=iT) \ 
VrR&idyT ~ dyS) 
iUi{So{y) + Q~iz=iT) 

Lemma 3.1. The operator Fi :V H is continuous. In particular, we have 

{Fi{u),u)=0. 
Proof. We have a constant ci > such that 
(9) IIV'llvi/|(D) < ci||g|lvi/,i(D) 

for any q E W2{D) which follows straightforwardly by regularity properties of a 
linear elliptic boundary problem. Note that W2 is a Sobolev space with respect to 
the third derivatives. Hence we get: 

||J(t,^)||l, < sup {\dyij{y,z)\ + \dMy,m^ 

{y,z)eD 



X I l^\dyT{y,z)\ + \dzT{y,z)\dD 
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The second factor on the right hand side is bounded by 

'J^\dyT(y,z)\^dDy + (^ljd,T{y,z)\'dDy < \\u\\v. 

On account of the Sobolev embedding Lemma, we have some positive constants 
C2, C3 such that 

sup i\dy^{y,z)\ + \d,ip{y,z)\) < C2\\Vi^\\wi{D) < cshWw^iD) < calMv 
Hence we have a positive constant C4 such that 

I|J(T»|U, <C4||U|12. 

for u G V. Similarly, we can treat the other terms containing J. 
We now show that 

(j(r,v),r)-o. 

For the other terms containing J we get a similar property. We obtain via integra- 
tion by parts 



dyTd.TpTdD- / d^TdytpTdD 

JD 

= - I dl^TipTdD+ I dlyT^TdD- [ dyT^pd^T dD + [ d^T^jdyTdD 

JD JD JD JD 

+ [ dyTi>T\lzldy - I d^TWllzldz = 

"'(0,1) "'(0,1) 

because ip is zero on the boundary dD. This relation is true for a set of sufficiently 
smooth functions V', T which are dense in wliD) x W2 (D). By^the continuity of Fi, 
as just shown in Lemma [3. II we can extend this property to W2{D) x W2{D). □ 

Lemma 3.2. The following estimate holds 

\\F2{u)\\v' < C5||w|k + C6. 

for some positive constants C5, Cg. 

Proof Let ( e W^{D). Sincc7;,=iT e H^O, 1) c 1), we have 7*^1 7^=1 T G 

(WliD))' and 

|(7:=i7.=iT,C)l = |(7.=iT,7.=iC)l < C7||T||^i||C||w^ 

for cy > 0, which immediately gives the first part of the above inequality. The other 
parts can be treated similarly. □ 

After this preparation, we are able to write our problem as a stochastic evolution 
equation. An introduction into the theory of stochastic differential equations one 
can be found in Zabczyk 33 . 

Let w be a noise on L2{D) with finite energy given by the covariance operator Q 
of the Wiener process w{t) which is defined on a probability space {n,J-,P). For 
the vector 

W = (u;, 0,0,0) 

we rewrite the coupled atmosphere-ocean system as a stochastic differential 
equation on V': 

nil 

(10) — + Au = F{u) + W, u{0) = uoeH, 
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where w is a white noise as the generalized temporal derivative of a Wiener pro- 
cess w with continuous trajectories on R and with values in ^2(0, 1). SufRcient for 
this regularity is that the trace of the covariance is finite with respect to the space 
i2(0, 1): tr^^Q < 00. In particular, we can choose the canonical probability space 
where the set of elementary events consists of the paths of w and the probability 
measure P is the Wiener measure with respect to covariance Q. 

In the following, we need a stationary Ornstein-Uhlenbeck process solving the linear 
stochastic equation on (0, 1) 

, , dz 

(11) d^ + ^i^-- 

where Ai = —dyy + {l + b(y)) is the linear operator with the homogeneous Neumann 
boundary condition at = and y = 1. 

Lemma 3.3. Suppose that the covariance Q has a finite trace : tr^^Q < 00. Then 
has a unique stationary solution generated by 

{t,Uj) ^ ziOtLo). 

Moreover, Z{lo) — (z(ti;), 0, 0, 0) is a random variable in V. 

For the proof we refer to Da Prato and Zabczyk j^Jj, Chapter 5, or Chueshov and 
Scheutzow 0. 

For our calculations it will be appropriate to transform H1Q|I into a differential 
equation without white noise but with random coefficients. We set 

(12) v:=u- Z 
Thus we obtain a random differential equation in V 

du 

(13) —+Av^Fi{v)+F2{v + Z{9tu;)), v{0) = vo e H. 
Equivalently, we can formulate the equation H13(l using test functions 

^^{v{t),C)+a{v{t)X) - {Fi{v{t))X) + {F2{v{t) + Z{0tu;)),O for aU C G V. 

We have obtained a differential equation without white noise but with random 
coefficients. Such a differential equation can be treated sample-wise for any sample 
u!. Hence it is simpler to consider (|13|) than to study the stochastic differential 
equation H10() directly. We are looking for solutions in 

V € C([0,r];i?) nL2(0,r; V), 

for all T > 0. If we can solve this equation then u :— v + Z defines a solution version 
of (|10|l . For the well posedness of the problem we now have the following result. 

Theorem 3.4. (Well- Posedness j For any time r > 0, there exists a unique 
solution of in C{[0,t]; H) D L2{0,t;V). In particular, the solution mapping 

R+ xnx H 3 {t,uj,vo) ^ v{t) e H 

is measurable in its arguments and the solution mapping 77 9 wq — > v{t) ^ H is 
continuous. 
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Proof. By the properties of A and Fi (see Lemma l3.1|l . the random differential 
equation is essentiaUy similar to the 2 dimensional Navier Stokes equation. 
Note that F2 is only an affine mapping. Hence we have existence and uniqueness 
and the above regularity assertions. □ 

On account of the transformation l|12() . we find that l|10|l also has a unique solution. 
Since the solution mapping 

IR+ X Q X H 3 {t,U!, vq) v{t, LU, vo) =: ip{t, lu, vq) G H 

is well defined, we can introduce a random dynamical system. On f2 we can define 
a shift operator 9t on the paths of the Wiener process that pushes our noise: 

w{-,9tuj) ^ w{- + t,uj) - w{t,uj) for t e R 

which is called the Wiener shift. Then {Ot}teR forms a flow which is ergodic for the 
probability measure P. The properties of the solution mapping cause the following 
relations 

(p(t + r, u!, u) = (p(t, 9rUj, ip{T, u!, u)) for t, T > 

(p{0, U!,u) — U 

for any oj € fl and u E H. This property is called the cocycle property of ip which 
is important to study the dynamics of random systems. It is a generalization of 
the semigroup property. The cocycle ip together with the flow forms a random 
dynamical system. 



4. DiSSIPATIVITY 

In this section we are going to show that the coupled atmosphere-ocean system 
is dissipative, in the sense that it has an absorbing (random) set. This definition 
has been used for deterministic systems . This means that the solution vector 
V is contained in a particular region of the phase space H after a sufficiently long 
time. Dissipativity will be very important for understanding the asymptotic dy- 
namics of the system. This dissipativity will give us estimate of the atmospheric 
temperature evolution under oceanic feedback. Dynamical properties that follow 
from this dissipativity will be considered in the next section. In particular, we will 
show that the coupled atmosphere-ocean system has a random attractor, has finite 
degree of freedom, and is ergodic under suitable conditions. 

We introduce the spaces 

H ^ L2{0,1) ^ L2{D) X L2{D) 
V = WUOA)^W^{D)xW^{D). 

We also choose a subset of dynamical variables of our system 

(14) i = {e,T,s), e^e-z. 
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To calculate the energy inequality for v, we apply the chain rule to HwHl/- We 
obtain by Lemma l3. II 

|||i;||| + 2||Vt)||i^+2||6^e||i^+2||e||i^ 

(15) ^-2ia,e)L,+2{Sa,e)L,-2{bSo,e)L,+2ibj,=iT,e)L, 

+ 2{^:^,So, T) + 2{^Ue, T) + 2(7:^iz(0t^), T) - 2{^U^.^,T, T) 
+ 2{jUF,S). 

Here and in the following we stress that < b{y) < 1. The expression Vw is defined 
by {dyQ,Vy^zT,Vy^zS). We now can estimate the terms on the right hand side. 
We have the following estimate for the second line of H15|l by the Cauchy-Schwarz 
inequality 



4a^ + Jll©llL+4||5J|i^ + i||e|,^^ 



+ nbisjl + IwbieWl^ + ^||6^||L^||7.=iT||i, + |||6^e| _ 
For the next line we obtain the estimate 
6||5ollL + i||7.=ir||L 

+ l|e||L + ll7.=ir||L + 6\mu;)\\l^ + i||7.=iT||L " 2||7.=iT||L. 

Now we can estimate the last line of H15|l . For any e > we can find an cs{e) > 
such that 

e\\Sr^r+cs{s)\\F\\l^. 

Here we also applied the trace theorem ||72=iS'|| i < C9||5||t4^i. Adding all terms 
containing ||7z=ir|||^^, we see that the sum is negative. 
Collecting all these estimates, we have 

' \\v\\l + 2\\vv\\l^ + |(i - ||6^|U^)||7.=iTiii2 + hewh 



dt' 

< 4a2 + 4||5„||i, + 2\\b^So\\l, + 6||^o|lL + 6||^(Mlli2 + cs{e)\\F\\l^. 

By using the Poincare inequality for S E W2{D), (Q, and choosing e small enough, 
we conclude that there is a positive dissipativity constanta such that 

d_ 
di' 

where a, cio is determined by physical data a^, ||S'a||L2i II'5'o||l2i II-^IU2 



(16) -\\v\\% + a{\\v\\l + W^^WV < + 6|k(et^)|lL 



By the Gronwall inequality, we finally conclude that 



(17) Ml < ||«(0)|||e-"* + ^ + 6e-"* f ||z(0,c^)||i^e"^d.. 

^ Jo 

We now show the dissipativity of v and v. Roughly speaking dissipativity means 
that all trajectories of the system move to a bounded set in the phase space. For a 
random system we have the following version of dissipativity. 
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Definition 4.1. A random set B = {B{Lj)}u)£n consisting of closed hounded sets 
B{uj) is called absorbing for a random dynamical system Lp if we have for any 
random set D — {_D(ti>)}ijgo, D{lo) G H bounded, such that t — > ^wpy^j^i^g^^-^ WdWh 
has a subexponential growth for t — > ±00 

ip{t,u,D{Lo)) C B{0tLo) fort>tQ{D,uj) 

v{t,e^t^,D{e.tio)) ^B{w) fort>to{D,Lu). 

B is called forward invariant if 

(p{t, uj, uq) g BiOtUj) if uq e B{uj) for t > 0. 

Although V is not a random dynamical system in the strong sense we can also show 
dissipativity in the sense of the above definition. 

Lemma 4.2. Let (p{t,U!,VQ) G H for vq E H be defined in U(J\) . Then the closed 
ball B{0, Ri{lu)) with radius 

Riiio) = 2 / e""(cio + 6\\z{erCj)\\lJdT 
J —00 

is forward invariant and absorbing. 

The proof of this lemma follows by integration of (|16|l . 

It remains to prove the dissipativity of the dynamical system Lp. To this end we 
obtain from the second equation of QJ: 

|ll'?llL+c?iPr||g|lL<PrRa2||V«||i, 



< 



PrRa^cio PrRa^ 2 PrRa^ d 2 

Q——\\z{Otu^)\\L. - ——yMh 



a a ^ a dt 

with sorne embedding constant cn in the Poincare inequality \\q\\L2 cii||Vg||L2 
for q G 14^2 (^)- Note that q satisfies homogeneous Dirichlet boundary conditions. 
Hence the variation of constants formula allows us to estimate: 



(19) 

" |f}(s)|||)e-^?iP^(*-^)ds. 



PrRa^ d 



a ds 

Now we apply the integration by parts 



ds 



e^-^'-*l|i)WliL + ll^(o)|lL 



Substituting this equation into ()19f) . we obtain 
(20) 

\\q{t,u;,u)\\l^ < ||z,o|||,e-=?iP'-*+ / (ilil^ cio + 6-^^ ||^(0,c.) ||i Je-^-P'-(*~^)ds 



* PrRa^ PrRa^ 
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Note that \\v(t)\\jj is bounded by 

which follows from H16() . To construct the radius of the absorbing set we have 
to replace uj by 6^t^. Suppose that t —> \\vo{9-t^)\\'jj growths not faster than 
subexponential. Then we have that 

lim ||i;o(w)||He-"* = 0, lim ||«o(0-t^)||l,e-"* = 0. 

t — >oc t — ^oo 

Hence we can conclude 



t 

and 



lim / £liil£_LL||„o(^)||2^e-"*e-=iiP'-(*-^)rfs = 







We also note that 



/■ ,PrRa Clo ^PrRa , ..o cfnRa Pr ^ ,^ _ 2 p,.f*_,-| , 

/■°PrRa2cio ,PrRa^, ,^ ,,,2 c^^Ra^Pr^ ,^ xNc^^Prsj 

The finiteness of this limit follows because the growth oit ^ Ri{9tuj) is subexpo- 
nential. Since all other terms are coupled with exponentially decreasing factors we 
have found: 



Lemma 4.3. Suppose that the assumptions of Lemma \4-<^ are satisfied. Then the 
random set {B{LLi)}i^^a given by closed balls B{0, R{uj)) in H with center zero and 
radius R{i-u) :— Ri{uj) + R2iuj) is an absorbing and forward invariant set denoted 
by B{uj) for the random dynamical system generated by fl!^) . 

For the applications in the next section we need that the elements which are con- 
tained in the absorbing set satisfy a particular regularity. To this end we introduce 
the function space 

■.= {ueH:\\u\\l:= < 00} 

where s G M. The operator A"^ is the s-th power of the positive and symmetric 
operator A. Note that these spaces are embedded in the Slobodeckij spaces H'^, s > 
0. The norm of these spaces is denoted by || • Wh"- This norm can be found in 
Egorov and Shubin 11 , Page 118. But we do not need this norm explicitly. We 
only mention that on the norm || • \\s of H'^ is equivalent to the norm of H'' for 
< s, see [H]. 

The reason to introduce these spaces is that the trace theorem can be formulated 
with respect to i?*; see Egorov and Shubin JT], Page 120: 

Theorem 4.4. Assume that a > ^. Then the trace mapping between two Sobolev/ Slobodeckij 
spaces 

is continuous. 
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This formula generalized Q because we can take for a = 1 and = W2{D) or V. 
Our goal is it to show that uj, D) is a bounded set in Ti^ for some s > 0. This 
property causes the complete continuity of the mapping w(l,a;, ■). We now derive 
a differential inequality for By the chain rule we have 

Note that for the embedding constant ci2,s between Ti" and V 




such that the left hand side is bounded if the initial conditions vq are contained in 
a bounded set in H. The second term in the above formula can be expressed as 
followed: 

t^{A^v,Aiv)H =2t{^v,A'v)H = -2t{Au,A''v)H + '2tiFi{v),A'v)H 
+ 2t{F2iv + Z{etij)),A'v)H. 

We have 

{Av,A%)h = \\A^+^v\\h = \\v\\1^,. 

If we apply some embedding theorems, see Temam Page 12 we have got for a 
Cl3 > 

Cl3 ll^l|mi + l ||'0||m2 + l llCllma C G -^ma 

where mi + m2 + ms > 1 and < to^ < 1. Here we use that D is of dimension 2. 
We then have for mi = 0, m2 — s < I and ma = 1 — s 

\{Fi{v),A'v)h\ < ci3\\v\\v\mi+,\\v\\i+s- 

11-0(^)111+5 is bounded by Ci||w(t)||ij similar to © and ||w(t)||L^(o,T;ff) < To 
ensure that the norm of ||i/i||m2+i is well defined we set ^ = (0, -0, 0, 0). Hence we 
have for any e > a constant ci4(e): 

(Fi(T;(t)),AMt))i^ <ci4(£)|k|lL(o,T;//)lkWlly + £||^WII?+s, 
where s is chosen sufficiently small. 

To obtain an estimate for the expression (^2(1; — Z),A^v) we only consider the 
expressions (7*^15*0,^*11) and {j*^ijz=iT, A^v) . The other term can be treated 
similarly. Suppose that S G ^2(0, 1). We interpret 7*5o as (0, 0, 0, 7*5o). Then we 
have for any e > a €15(2), Ci6(£) > such that for < s < i 

{^USo,A%)<s\\Ai+^v\\j, + ci,ie)\\A-i+ir^^,So\\j, 

<e\\v\\l+, + cieie)\\So\\l,_i < e\\v\\l+, + cie{e)\\So\\l^^oA)- 

Since by the trace theorem we have 7^=1 : ^^^+"(0, 1) iJ^^^", 1 — s > i. For 
a < the space H°' denotes the dual space of H^". 



14 



JINQIAO DUAN, HONGJUN GAO, AND BJORN SCHMALFUSS 



For {'^*^ijz=iT, A^v) we have for positive constants C17-C19, a sufficiently small 
e > 0, C2o(£) and a sufficiently small e' > 

(7:=i7z=iT, A'v) Ml.^iT, 7.=i^^«)l, 

< ||7.=ir|UJ|7.=i^^^^lU. 

< Ci7\\v\\v\hz=lA'v\\L, 

< Cis\\v\\v\\Ai+''A'v\\H < Cig\\v\\v\\Ai + iv\\H < Ci9||w||y 

<C2o(e)||«fy+e||«||?+,. 

For < s < i we have A^v G H^"'^" for w e F such that j^^iA^v G ^2(0,1). 
Collecting all these estimates we obtain that t\\v{t,uj,VQ)\\ff is bounded for t <T < 
00 if Wo is contained in a bounded set. This allows us to write down the main 
assertion with respect to the dissipativity of this section. 

Theorem 4.5. For the random dynamical system generated by there exists a 

compact random set B = {B{uj)}^^fi which satisfies Definition \4-l\ 

We define 

(21) B{cj) = ^(1, e^iuj, 5(0, R{0-iuj))) c n", 0<s<i 

In particular, Ti" is compactly embedded in H. 

5. Random dynamical behavior 

In this section we will apply the dissipativity result of the last section to analyse 
the dynamical behavior of the coupled atmosphere-ocean system However, it 
will be enough to analyse the transformed random dynamical system generated by 
(|l()|l . By the transformation (|12(l we can take over all these qualitative properties 
to the system (jlOII . 

We will consider following dynamical behavior: random climatic attractors, finite 
degrees of freedom, atmospheric temperature evolution under oceanic feedback, and 
ergodicity. 

We first consider random climatic attractors. We recall the following basic concept; 
see, for instance, Flandoli and SchmalfuB 12 . 

Definition 5.1. Let tp be a random dynamical dynamical system. A random set 
A = {^(w)}^gsi consisting of compact nonempty sets A{uj) is called random global 
attractor if for any random bounded set D we have for the limit in probability 

(P) lim distHMt,uj,D{Lj)),A{etUj)) = 

t — *oo 

and 

ip{t,u;,A{uj)) = A{9tu;) 

any t > and lu d fl. 

The essential long-time behavior of a random system is captured by a random 
attractor. In the last section we showed that the dynamical system (p generated by 
(|l()|l is dissipative which means that there exists a random set B satisfying p8(l . In 
addition, this set is compact. We now recall and adapt the following theorem from 
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Theorem 5.2. Let ip be a random dynamical dynamical system on the state space H 
which is a separable Banach space such that x — > Lp{t, uj, x) is continuous. Suppose 
that B is a set ensuring the dissipativity given in definition \4.1\ In addition, B has 
a subexponential growth (see Definition \4-l}) and is regular (compact). Then the 
dynamical system ip has a random attractor. 

This theorem can be appUed to our random dynamical system ip generated by the 
stochastic differential equation HlOfl . Indeed, all the assumptions are satisfied. The 
set B is defined in Theorem 14.51 Its subexponential growth follows from B{ij) C 
B{Q,R{i^)) where the radius R{i^) has been introduced in the last section. Note 
that ifi is a continuous random dynamical system; see Theorem 13.41 Thus ip has 
a random attractor. By the transformation (|12|l . this is also true for the original 
coupled atmosphere-ocean system. 

Corollary 5.3. (Tlandom Attractor j The coupled atmosphere- ocean system 
has a random attractor. 



Dissipative systems often have finite degrees of freedom. This is reflected by the 
fact that the Hausdorff-dimcnsion of the attractor is finite. This fact can be applied 
to fluid dynamical systems; see for instance Temam PO], Page 403 ff. A similar the- 
ory has been developed for random dynamical systems; see, for example, [71 IHl I23j . 

However, we will follow another approach to show that the random climatic at- 
tractor of has only finitely many degrees of freedom. Namely we will use the 
technique of determining functionals. This technique has been introduced for de- 
terministic systems by Foias and Prodi ^31 , and Ladyzhenskya . See ^| ^ for 
more recent work. Roughly speaking, a set of determining functionals is a set of 
functionals (for instance, Fourier modes) such that if it is known that a dynamical 
system has an asymptotic stable behavior only with respect to these finitely many 
modes, then the complete system has an asymptotic stable behavior. For random 
dynamical systems we can investigate determining functionals in probability as in 
Chueshov et al. 

Definition 5.4. We call a set C — {Ij, j = 1, ■ • • , N} of linear continuous and 
linearly independent functionals on a space X continuously embedded in H (for 
instance X = or V ) asymptotically determining in probability if 

rt+l 

(P) lim / max |L-((^(r, w, ui) — (^(r, w, U2))pdT = 

for two initial conditions Vi, V2 ^ H implies 

(P) lim \\ip{t,uj,vi) ~ ip{t,uj,V2)\\H = 

t — )-oo 

Often the elements of set C can be chosen as the projections with respect to the 
Fourier expansion of the solution. Since our domain D is an rectangle, we can 
calculate the Fourier expansion more explicitly. 

In the following we need an additive embedding inequality based on qualitative 
difference of the spaces H and X for some set of functionals C 



(22) 



\\u\\h < Ccma.x\k{u)\+ec\\u\\x, Cc > 
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The constant ££ > describes a fundamental difference of the spaces X and H. 
How (j221) works is described by an motivating example in Chueshov et al. [H]. 
We recall and adapt a result from [S]. 

Theorem 5.5. Let C = {Ij : j — 1,...,A^} be a set of linear continuous and 
linearly independent junctionals on X. We assume that we have an absorbing 
and forward invariant set B in X such that for supj,^^^^) ||u||^ the expression 
t — > s\vp^i^B(etuj) Ikllx locally integrable and subexponentially growing. Suppose 
there exist a constant C21 > and a measurable function I > such that for 
wi, V2 V we have for F{u!, v) — Fi{v) + F2{v + Z{u!)) 

{-A{vi - V2) + F{uj, vi) ~ F{uj, V2), vi - V2) 

< -C2l\\vi - V2\\v + l{vi,V2,Uj)\\vi - V2\\h- 

Assume that 

— E< sup / l{if{t,u},vi),ip{t,u},V2),9tuj)dt\<C2ie]-^^ 

ITT- [vuV2eB{ui)Jo J 

for some m > 0. Then C is a set of asymptotically determining functionals in 
probability for random dynamical system if. 

We set X — H^,s £ (0, j). In the last section we have shown that the set B, 
consisting of bounded sets, is forward invariant. The function / appearing in the 
formulation of Theorem IS . 51 expresses the local Lipschitz continuity of our nonlinear 
operator F{v). The essential part that determines F is the Jacobian operator J. A 
method on how to estimate the function / is in jSj . We also note the local Lipschitz 
constant can be estimated in terms of the TY'^-norm. To estimate F2 we have to 
apply the techniques introduced in Section 4. Note that F2 is linear. 
We can apply Theorem 15. 51 to the random dynamical system generated by pU|l and 
get the following result. 

Theorem 5.6. The random dynamical system generated by has finitely many 
degrees of freedom. More precisely, there exists a set of linearly independent con- 
tinuous functionals (on TL'^ ) which is asymptotically determining in probability. 

Because of the transformation (|12|l . this result is also true for the coupled atmosphere- 
ocean system: 

Corollary 5.7. (T'inite Degrees of Freedomj The coupled atmosphere- ocean 
system 0) has finitely many degrees of freedom, in the sense of having a finite set 
of linearly independent continuous functionals which is asymptotically determining 
in probability. 

The number of the elements in C can be estimated explicitly as in [2| . 

Now we consider random fixed point and ergodicity. We can do a small modification 
of (Q. This modification is given when we replace AT by j/AT and AS* by v/S.S 
where > is viscosity. Under particular assumptions about physical data in , 
we can show that the behavior of our dynamical system is laminar. For a stochastic 
system, this means that after a relatively short time, all trajectories starting from 
different initial states show almost the same dynamical behavior. This can be seen 
easily if ^o, a, Sa, F are zero, there is no noise and v is large. We will show that a 
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laminar behavior also appears when So, a, Sa, F are small in some sense. 
Mathematically speaking, laminar behavior means that a random dynamical system 
has a unique exponentially attracting random fixed point. 

Definition 5.8. A random variable v* : CI ^ H is defined to be a random fixed 
point for a random dynamical system if 

ipit,iU,V*{Lu))^V*{etu) 

for t > and lo Q. A random fixed point v* is called exponentially attracting if 

lim \\ipit,iu,x)-v*{etu;)\\H = 

t — >oo 

for any x Cz H and uj E Q. 

Sufficient conditions for the existence of random fixed points are given in Schmalfuf5 
|24| . We here formulate a simpler version of this theorem and it is appropriate for 
our system here. 

Theorem 5.9. (Handom Fixed Point Theorem^ Let (p be a random dynamical 
system and suppose that B is a forward invariant complete set. In addition, B has 
a subexponential growth, see Definition \4-l\ Suppose that the following contraction 
conditions holds: 

(23) sup < ^(^) 

where the expectation o/logfc denoted by Elogfc < 0. Then ip has a unique random 
fixed point in B which is exponentially attracting. 

This theorem can be considered as a random version of the Banach fixed point 
theorem. The contraction condition is formulated in the mean for the right hand 
side of 



Theorem 5.10. Assume that the physical data \a\, WSaWb^i II*5'o||l2j II-^IU2 '^"'^ 
the trace of the covariance for the noise tinQ are sufficiently small, and that the 
viscosity v is sufficiently large. Then the random dynamical system generated by 
^1(J\) has a unique random fixed point in B. 

Here we only give a short sketch of the proof. Let us suppose for a while that B is 
given by the ball B{0,R) introduced in Lemma r4.3l Suppose that the data in the 
assumption of the lemma are small and v is large. Then it follows that ER is also 
small. To calculate the contraction condition we have to calculate uj, vi{lu)) — 
(^(1, Lu, V2{Ltj))\\'j^ for arbitrary random variables vi, V2 G B. By Lemma l3. II we have 
that 

{J{qi,tjji) - J{qi,ipi),qi - gs) < C22\\qi - q2\\w^ + C23||gi||^i(r,)||gi - 92||i2 

where the constant C23 can be chosen sufficiently small if v is large. On account of 
the fact that the other expressions allow similar estimates and that F2 is linear we 
obtain: 

^\\ip{t,UJ,Vl{uj)) - ip{t,UJ,V2{uj))\\H 

< (-a' + C23\\(p{t, UJ, V2(uj))\\v)\\(p{t, W, Vi{uj)) - ip{t, UJ, V2{uj))\\h 
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for some positive a' depending on v. From this inequality and the Gronwall lemma 
it follows that the contraction condition H23(l is satisfied if 

E sup C23 /" \\'f{t,uj,V2)\\vdt < a' . 

U2£B{uj) Jo 

But by the energy inequality this property is satisfied if the Ei? and E||z|jy is suf- 
ficiently small which follows from the assumptions. 

Let now B be the random set defined in |2U- Since the set B introduced in pifl 
is absorbing any state the fixed point v* is contained in this B. In addition v* 
attracts any state from H and not only states from B. 

Corollary 5.11. (Unique Random Fixed Point j Assume that the physical 
data\a\, \\Sa\\L2^ II'5'o||l2i II-P'I|l2 and the trace of the covariance for the noise trnQ 
are sufficiently small, and that the viscosity v is sufficiently large. Then, through 
the transformation f Jjj)). the original coupled atmosphere- ocean system 0) has a 
unique exponentially attracting random fixed point u*{lo) = v*{lj) + Z{lo), where 
u = {Q,q,T,S). 

The uniqueness of this random fixed point implies ergodicity. We will comment on 
this issue at the end of this section. 

By the well-posedness Theorem 13.41 we know that the stochastic differential equa- 
tion I|1U|) for the coupled atmosphere-ocean system has a unique solution. The 
solution is a Markov process. We can define the associated Markov operators T(i) 
for t > 0, as discussed in |21||111. Moreover, {T{t)}t>o forms a semigroup. 
Let be the set of probability distributions /i with finite energy, i.e., 

< OO. 

Then the distribution of the solution u{t) (at time t) of the stochastic differential 
equation H1U|) is given by 

where the distribution fiQ of the initial data is contained in il/^. 

We note that the expectation of the solution ||w(i)|||^ can be expressed in terms of 

this distribution T(t)fiQ: 

nHt)\\l= I \\ufHdT{t)^io. 

J H 

We can derive the following energy inequality in the mean, using our earlier esti- 
mates in and ISOJ: 

Theorem 5.12. The dynamical quantity u — (0, q^ T, S) of the coupled atmospheric- 
ocean system satisfy the estimate 

E||w(t)|||,+aE / \\u{T)\\ldT <¥.\\u4]j + tc2A + tirL2Q, 
Jo 

where the positive constants C24 and a depend on physical data F{y), a, Sa{y), 
So{y), Pr and Ra. 

By the Gronwall inequality, we further obtain the following result about the as- 
ymptotic mean-square estimate for the coupled atmosphere-ocean system. 
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Corollary 5.13. (^Atmospheric Temperature Evolution under Oceanic Feed- 
back^ For the expectation of the dynamical quantity u = (0, T, S) of the coupled 
atmospheric- ocean system, we have the asymptotic estimate 

irL2Q 



limsupE||u(t)||^ = limsup / ||u||^dr(t)^o < 



C25 



if the initial distribution fio of the random initial condition mo(^) is contained in 
M'^ . Here C25 > also depends on physical data. In particular, we have asymp- 
totic mean-square estimate for the atmospheric temperature evolution under oceanic 
feedback 

(24) limsupElieill, < + 

t-^oo C25 

Thus the atmospheric temperature Q{y,t), as modeled by the coupled atmosphere- 
ocean system is bounded asymptotically in mean-square norm in terms of phys- 
ical quantities such as the freshwater flux F{y), the trace of the covariance operator 
of the external noise, the earth's longwave radiative cooling coefficient a, and the 
empirical functions Sa{y) and So{y) representing the latitudinal dependence of the 
shortwave solar radiation, as well as the Prandtl number Pr and the Rayleigh num- 
ber Ra for oceanic fluids. 

By the estimates of Theorem 15.121 we are able to use the well known Krylov- 
Bogolyubov procedure to conclude the existence of invariant measures of the Markov 
semigroup. 

Corollary 5.14. The semigroup of Markov operators {T(t)}t>o possesses an in- 
variant distribution in : 

T{t)^i = for t>0. 

In fact, the limit points of 

J J T{T)nodT\ 

for t 00 are invariant distributions. The existence of such limit points follows 
from the estimate in Theorem 15. 121 

In some situations, the invariant measure may be unique. For example, the unique 
random fixed point in CoroUarv 15.111 is defined by a random variable u*{uj) = 
v*{uj) + z{ll!). This random variable corresponds to a unique invariant measure 
of the Markov semigroup. More specifically, this unique invariant measure is the 
expectation of the Dirac measure with the random variable as the random mass 
point 

Because the uniqueness of invariant measure implies ergodicity |22) . we conclude 
that the coupled atmosphere-ocean model (Q) is ergodic under the suitable con- 
ditions in Corollary 15.111 for physical data and random noise. We reformulate 
CoroUarv 15. Ill as the following ergodicity principle. 

Theorem 5.15. (^Irgodicityj Assume that the physical data \a\, \\Sa\\L2J II'S'o||l2: 
||-F||l2 cmd the trace of the covariance for the noise trnQ o-fs sufficiently small, and 
that the viscosity v is sufficiently large. Then the coupled atmosphere- ocean system 
iPP is ergodic, namely, for any observable of the coupled atmosphere-ocean flows, 
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its time average approximates the statistical ensemble average, as long as the time 
interval is sufficiently long. 

6. Summary 

We have investigated the dynamical behavior of a coupled atmosphere-ocean model. 
First, we have shown that the asymptotic dynamics of the coupled atmosphere- 
ocean model is described by a random climatic attractor fCorollarv l5.3|l . Second, 
we have estimated the atmospheric temperature evolution under oceanic feedback, 
in terms of the freshwater flux, heat flux and the external fluctuation at the air-sea 
interface, as well as the earth's longwave radiation coefficient and the shortwave 
solar radiation profile fCoroUarv 15. f 3(1 . Third, we have demonstrated that this 
system has finite degree of freedom by presenting a finite set of determining func- 
tional in probability fCoroUarv I5.7|l . Finally, we have proved that the coupled 
atmosphere-ocean model is ergodic under suitable conditions for physical parame- 
ters and randomness, and thus for any observable of the coupled atmosphere-ocean 
fiows, its time average approximates the statistical ensemble average, as long as the 
time interval is sufficiently long fTheorem l5.15(l . 
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